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1 Introduction 



Pontecorvo's insightful idea of neutrino oscillations [[[[] has met phenomenal success over the 
past several decades in explaining the observed deficit of neutrinos from a wide variety of 
astrophysical and terrestrial sources. It is now well established that the different neutrino 
flavors oscillate among themselves while conserving the lepton number. The coherent phase 
picked up by the evolving states that shows up in transition amplitudes leads to the phe- 
nomena of oscillations. The fact that there is a geometric interpretation of the standard 
neutrino oscillation formalism was recently shown for the specific case of CP conservation 
(CPC) and two flavors. 

The topological 1 phase obtained in Ref. was robust and did not depend upon the 
nature of the evolution. We pointed out that as long as the singular point in ray space 
was encircled, the phase of ir remained irrespective of the evolution being in vacuum or in 
medium. The ir anholonomy was first noticed in Refs. and [|J as a phenomenon associated 
with real and symmetric Hamiltonians. The presence of this topological phase was ascribed 
to the structure of leptonic mixing matrix. This implied that there was no unexpected 
surprise due to this topological phase and that the standard quantum mechanical treatment 
takes into account Pancharatnam's phases correctly (in fact is a realization of the same). 

The purpose of the present work is to extend our analysis || to the CP violating (CPV) 
case while restricting ourselves to two neutrino flavors. We will explicitly find the conditions 
under which it is possible to destroy the robustness of the topological phase by the inclusion 
of CPV terms. We consider the general CPV form of neutrino Hamiltonian for the case of 
two neutrino flavors and derive an expression for the oscillation probability for an arbitrary 
medium with in the adiabatic approximation. Then we show that the cross terms in the 
probability have a purely geometric interpretation by employing Pancharatnam's prescrip- 
tion of cyclic quantum collapses. By doing so we establish a crucial geometric difference 
between the CPC and the CPV case. Having established our result, we briefly discuss our 
result in two possible physical scenarios where CPV can appear in the Hamiltonian and 
leave its imprint at the level of detection probability. 

The information on neutrino masses and mixings gleaned from most oscillation ex- 
periments can be conveniently analysed considering two neutrino situation or quasi-two- 
neutrino-situation (for example, the one mass scale dominant approximation) ||. So, we 
focus only on two flavor case in the present work which provides a clean visualization of 
physical effects via rotations and reflections on the Poincare sphere. 

This paper is organized as follows. In Sec. |2|, we describe the general Hamiltonian for 
two flavor neutrinos, including the CPV term. We also discuss possible sources of such 
a CPV term. Then in Sec. [3], we obtain the neutrino oscillation formulae in presence of 

1 By our definition, the topological phase refers to phase factors that are insensitive to small changes in 
the circuit (and are invariant under deformations of circuit), while geometric phases are sensitive to such 
changes. We use this distinction between the two terms throughout this article (unless otherwise specified). 
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CPV term and establish connection of the cross terms appearing in the probability with 
the Pancharatnam's geometric phase. We explicitly demonstrate (pictorially) the geometric 
differences between the case of constant density and adiabatically changing matter back- 
ground, both for CPC and CPV situations which is the central result of this article. In this 
context, we also comment on the action of CP and its connection with Wigner's theorem in 
Sec. f|. We end with concluding remarks in Sec. [5]. 



2 Hamiltonian for two flavor neutrinos 

It is well-known that the two flavor neutrino system is equivalent to a two state quantum 
system 0, ^| in the ultra-relativistic limit (for equal fixed momenta of the two neutrinos). 
Upon linearizing the Dirac equation in the ultra-relativistic limit, the propagation of neutri- 
nos through vacuum or matter can be described by an effective Schrodinger-like equation. 

Let us examine the general hermitian 2 form of Hamiltonian for any two state system , 

H = f ■ a + rob 

( z x — iy\ 
\x + iy — z J 

^ (— cos-# sin?? e~ lip \ , . 

~ ^sintf e** cos$ J ' ^ ' 

where, x, y and z are the three independent parameters appearing in the Hamiltonian. 
The term proportional to Identity r leads to non-zero trace and adds an overall phase 
factor, e~ l f r ° dt to the evolving state. But the states have the freedom of redefinition upto 
a phase so this extra phase does not affect oscillations and we can safely omit this term 
and deal with traceless Hamiltonian (last row of Eq. ([!])) spanned by the Pauli matrices. 
$ and ip are the polar angles and r is a unit vector. The eigenvalues of Eq. ([!]) are A± = 
± J Q * a/x 2 + y 2 + z 2 dt' = ±1 and the normalized eigenstates are : 

, Q . /cos(tf/2)e-^\ , , . > /-sin(tf/2)e-^\ 

l*'*. + >=( ) and l^^->=( cosW2) J' ^ 

The ray space in this case is a two-dimensional sphere S 2 which is known as Poincare sphere 
(also referred to as the Bloch sphere in quantum mechanics). Pictorially, oscillations are 
unitary rotations about an axis on the Poincare sphere while CP transformation can be 
represented as reflection in the x — z plane on the Poincare sphere. The terms appearing in 
the H due to vacuum, normal matter (with standard interactions (SI) and non-standard in- 
teractions (NSI)) and due to neutrino backgrounds with SI encountered by the neutrinos as 

2 neglecting absorption 
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Medium 


X 


y 


z 


Vacuum 


(w/2) sin -d 





-(w/2) cost? 


Normal 
matter+SI 


(w/2) sintf 





-(uj/2) cos-d + V c /2 


Normal 
mattcr+NSI 


Re[(w/2)sin*+e e «/2] 


Im \{ui/2) sintf + e e „/2] 


-(u/2) cos-d + V C /2 + {e ee - e yy )/2 


Neutrino 
backgrounds+SI 


Re[(w/2)sini? + B e!/ /2] 


Im[(w/2) sintf + B ey /2] 


~{u)/2) cos-d + V c /2 + B/2 



Table 1 : The three independent elements of H in different kinds of media - vacuum, normal matter with 
SI, ordinary matter with NSI and mixed state neutrino backgrounds (see text). 



they propagate between the source and the detector are listed in Table |ip. In general, oscil- 
lations in vacuum or any medium can be viewed as the phenomenon of elliptic birefringence 
in optics as explained in 0. 

Let us comment upon the CP properties of different media mentioned in Table || and 
appearance of CPV term in the Hamiltonian. The action of CP can be understood in terms 
of H — > H T . Since only a y changes sign under CP, CPV implies that y ^ 0. Vacuum is CP 
symmetric which implies y = 0. The other two terms (x and z) depend on mixing angle 
in vacuum d/2 and the vacuum oscillation frequency u = 5m 2 /2p where 5m 2 is the mass- 
squared difference and p ~ E is the neutrino momentum or energy in the ultra-relativistic 
limit. Normal matter is CP asymmetric since neutrinos and antineutrinos interact with 
matter in different ways (Mikheyev-Smirnov-Wolfenstein || Q effect can enhance oscillation 
in the neutrino channel and suppress in the antineutrino channel) but due to the absence 
of flavor changing neutral currents in the Standard Model (SM) there are no off diagonal 
terms and the Hamiltonian still preserves CP (y = 0) in the absence of intrinsic CPV phases. 
An important consequence of this is that if vacuum mixing and neutrino masses are zero, 
then neutrinos can not change flavor. In Table [1], the second row corresponds to this case. 
Vc = V2Gpn e is the charged current potential due to coherent forward scattering of v e 
with electrons in matter. Gp is the Fermi coupling constant and n e is the electron number 
density. 

The minimal theoretical scenario needed to describe neutrino flavor oscillations is the 
requirement of neutrino masses, for which there is concrete evidence from various oscilla- 
tion experiments. It is well-known that the SM by itself predicts massless neutrinos. The 
simplest way is to add dimension five non-renormalizable terms consistent with symmetries 
and particle content of the SM which leads to desired Majorana masses for the left-handed 
neutrinos. But, the neutrino interactions involving the light fields are still assumed to be 
described by weak interactions within the SM and this is minimal SM extension for ac- 
comodating neutrino mass. However, once we invoke new physics in order to explain the 

3 The angle d is connected to the mixing angle O used in the standard neutrino oscillation literature by 
t9 = 26. 
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non-zero neutrino masses, it seems rather unnatural to leave out the NSI which allow for 



flavor changing interactions as well as are new sources of CP violation which can affect 



production, detection and propagation of neutrinos [10], |TT[. Some of the early attempts 
discussing new sources of lepton flavor violation (for instance, R-parity violating supersym- 
metry) were geared towards providing an alternate explanation for the observed deficit of 
neutrinos coming from the Sun in the limiting case of zero neutrino mass and absence of 
vacuum mixing [12, 13, 14]. In recent years, the emphasis has shifted towards understanding 



the interplay between SI and NSI and whether future oscillation experiments can test these 
NSI apart from determining the oscillation parameters precisely. This has led to an upsurge 



in research activity in this direction [[15 , 16|, |17| , hq , [L9J f2(| |2T], 0, [23] . Non-standard physics 



which allows for flavor changing interactions may also lead to additional CPV phases which 
appear as y ^ in Eq. ([!]). 

The first possibility to have y ^ is the case of NSI in propagation (which affects the 
coherent forward scattering of neutrinos with background matter) 4 which induces new CPV 
phases. Such extrinsic 5 CPV phases can mimic the signature of intrinsic CPV phases thereby 
affecting the determination of intrinsic phase appearing in the leptonic mixing matrix. In 
the three flavor case, the problem becomes very complicated and one needs to carefully 
disentangle the extrinsic CPV from the intrinsic one by looking at the CPV observables [|9[ 



Here in the two flavor case, there are no intrinsic phases |24J but extrinsic CPV phases can 
actually appear at the level of Hamiltonian. If we consider only those operators that arise 
at a scale much lower than lepton number violating scale then the relevant interaction is 

£ = J2 ^^VaLl^vpL UiphljL + e f Jf R ^f R ) , (3) 

where / = e,p,n and a, (3 = e,fi,r. Since the scale at which this interaction arises is 
supposed to be not too far from the electroweak scale A 

ew ~ G p , its coupling may 
be parameterized by e Q( gGV, where e ~ (A^/A^p) 2 . The terms appearing in the flavor 
Hamiltonian in presence of ordinary matter and NSI are given in third row of Table |l[ At 
the level of underlying Lagrangian describing NSI, the NSI coupling of the neutrino can be 
to e,u,d or e,p,n. But from phenomenological point of view, only the sum (incoherent) 
of all these individual contributions is relevant. The effect of coherent forward scattering 
induced by such interaction terms in Lagrangian on the propagation of neutrinos in ordinary 
neutral unpolarized medium is governed by the parameters, 



n 



f=e,u,d 

where nj is the density of fermion / in medium crossed by the neutrino. Also, e* = e^ L -\-e' R . 



4 We focus only on NSI during propagation of neutrinos through matter. However it should be noted 
that NSI can also affect production and detection processes as discussed in Refs. [|5[ [l6[ ^ g3[ . 
5 The term extrinsic refers to matter-induced CPV phases. 
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Note that a CPV term in the Hamiltonian can have its origin in different sources, 
for example by inclusion of non-standard neutrino-matter interactions during propagation 
as described above or via CPV neutrino background which is encountered by neutrinos 
emanating from a dense supernova core. In the latter case, if the background neutrinos are 
in a mixed 6 state given by \i*b) — 7e| v e ) + ly \ v y } , CPV can arise in the two flavor neutrino 



Hamiltonian without the need to invoke any new physics [gg, ^6|. This is due to neutrino 
self interactions [27, 28, p9[] . The pure state flavor density matrix for a given momentum 
mode p is 

ft. = I = y$) ■ ( 5 ) 

Background neutrinos will be in mixed state naturally if they also undergo oscillations. The 
terms of fourth row in Table p] are given by 

B = V2G F J d 3 q(l - cos# pq )[(£ q - £? q ) ee - (g q - Q^ yy ] , 

B ey = V2G F J rf 3 q(l - cos6> pq )[(£ q - Q^) e y] , 

B ye = \PlGp J t/ 3 q(l COS$ pq )[(£> q f? q )ye] i (6) 

where # pq is the angle between the direction of propagating neutrino with momentum p 
and the direction of other neutrinos in the ensemble with momentum label q. The quantity 
Qqey^Qqeu) stands for the matrix element of the density matrix operator, ( v e | £> q | v y ) (( v e \ £> q | v y 
Hence for mixed state background neutrinos, all the three terms of the Hamiltonian x, y and 
z can be non-zero in contrast to the case of vacuum or ordinary matter with SI. Thus we 
note that the Hamiltonian with CPV neutrino background requires all the three parameters 
for its complete description. The background potential matrix is no longer diagonal and 
therefore even massless neutrinos can oscillate. 

We should emphasize that within the SM, the only situation where one can get flavor 
off-diagonal terms in the Hamiltonian without any new physics is the case of dense neutrino 
backgrounds. Here off-diagonal terms arise not due to any new flavor changing interactions 
but due to background neutrinos being in a mixed flavor state. If the background neutrinos 
were in pure flavor state then the off-diagonal terms in the potential matrix would disappear 
and one ends up with a CPC Hamiltonian as in the case of ordinary matter when only SI 
are taken into account. Hence such a mixed state neutrino background provides us with a 
minimal scenario for CPV effects to play a role. In addition, invoking NSI in the case of 
dense neutrino backgrounds leads to extra source of CPV phases in the Hamiltonian |EU 
over and above those induced by mixed state neutrino background. 

The question then arises if there is any visible consequence of extrinsic CPV phase or 
y ^ at the level of probability in the two flavor case whatever its source might be. It 



3 Not an impure state in the conventional sense. 
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is well-known that in the case of CPV medium with constant density, this term is of no 
consequence. Below we will first establish this fact by showing that the extra phases in the 



mixing matrix can be gauged away |24|. By using Pancharatnam's prescription of cyclic 
quantum collapses we show that the geometric phase is always restricted to be topological 
(two pairs of orthogonal states can always be made to lie on a great circle) || . Further we will 
discuss the case of varying density CPV medium where it is possible to have an observable 
effect of CPV phases at the probability level. This is connected intimately to the fact that 
the Pancharatnam's phase can become geometric (different from it, 0) in such a situation. 
We give a pure geometric view of this effect by using Pancharatnam's prescription [R(|. 



3 CPV phases and the two flavor oscillation probabil- 
ity 

Solving the full neutrino evolution equation in matter exactly is formidable analytically 
because of the fact that density in general is not constant but varies in space (or time in 
natural units). Usually one assumes adiabaticity : the variation of density is small over one 



oscillation length. This is the quasi-static approximation pT|. 

If we start with the most general form of mixing matrix which includes the CPV phases, 
we can argue that only one angle is sufficient to parameterize this matrix and the extra phase 



can always be gauged away by appropriate redefinition of the flavor states [gj, |32], [33[. This 
can be shown explicitly as follows 



cos(tf/2)e- sin(0/2)e* 
" I - sin(i?/2)e*T cos(^/2)e J (- a+ ' 3+ ^ 



(7) 



If we define 



p » - y e -*y J ' Pi - y e i(-a+p) J » ( 8 ) 

then under a transformation U — > P U VP^, we can eliminate the three phases from the mixing 
matrix. The mass eigenstates are defined up to a phase and therefore we can redefine them 
u l,r Pv v l,r (Note that u L is denoted by | tfj, (ft, ± ) where 2 = 1,2 runs over the number 
of generations) and similarly for l L R — > PiIl,r- With the above transformation, the mass 
matrices remain real and unchanged. This leads to no observable CPV (Dirac type) in 
the two flavor case. But, a Majorana type CPV phase can survive even in the two flavor 
B2], B3, B3, B5fl. Note that in certain extended models, such as left right symmetric 



case 



models and supersymmetric models, it is possible to have CPV only with two flavors. The 
CPV characteristic of Majorana phases can be probed in effects connected with the neutrino 
mass term such as neutrinoless double beta decay MM . But the main point is that neutrino 
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oscillation experiments are insensitive to Dirac or Majorana nature of neutrinos |32|, [35 
So, it suffices to take the mixing matrix to be real, 



fcos(0/2) sin(tf/2)\ (q) 
U V" sin W 2 ) cos($/2)J ' 1 ' 

without any loss of generality. This leads to no observable significance of intrinsic CPV 
phases as far as two flavor neutrino oscillation formalism is concerned PBJ. Note that at 



any instant, we can always get rid of the intrinsic CPV phases in the mixing matrix. 

However in case of an adiabatic evolution, is it still possible to define a gauge transfor- 
mation which can gauge away the dependence of CPV phases. Let us ask if the transition 
probability for the two flavor neutrinos is at all sensitive to the CPV phases in presence of 
adiabatic evolution. In what follows, we consider neutrinos propagating through an arbi- 
trary medium and obtain the expression for transition probability and analyse the effect of 
CPV phase in the Hamiltonian in Eq. ([I]). 

To keep the discussion general, we start with a neutrino created as a flavor state | v a ) 
and detected as \vp). The state \v a ) is 

| v a ) = v aj+ \'&i,ip 1 , + ) + u a -\ &i,<px, -> , (10) 

where | ipi, ± ) denote orthogonal pair of energy eigenstates of the Hamiltonian H($i, <pi) 
(Eq. (JJ)). The two orthogonal mass states | $i, <pi, ± ) evolve into | $2, y?2, ± ) as a result of 



adiabatic evolution [EH] 



0i,pi,±) e~ lV± \d 2 ^2,±) , (11) 



with V± = ± J Q \/ x 2 + y 2 + z 2 dt' + J Q r dt' as the dynamical phases obtained from Eq. ([I]) 
where x, y, z for specific medium can be read off from Table [1]. 

This is the zeroth order adiabatic approximation. However the two states can also 
pick up a geometric contribution as they traverse the region between the source and the 
detector 0. In general the geometric component picked up by the individual eigenstates 
under adiabatic evolution vanishes on account of evolution being along the geodesic curve 7 . 
But, a net geometric phase can appear if we can perform interference experiment between 
these energy eigenstates in the energy space. Neutrinos are produced and detected as flavor 
states which are superpositions of mass states, so this condition of realising interference 
in energy space is automatically achieved [0]. The appearance of geometric phase can be 
understood from the description of vertical lift |36] which accounts for the net non-zero 
geometric phase during any evolution (closed or open). Note that open loops can also be 
closed by connecting them by the shorter geodesic. The geometric phase is determined by 
the path of the state vector on the Poincare sphere. 



7 Except when the evolution along the geodesic curve happens to cross the antipodal point and the 
geometric phase then becomes tt. 
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The states | tpi, ± ) and | $2, <^2, ± ) are connected by parallel transport rule on the 
Poincare sphere. The two time evolved states e~ lT>± \ $2, <P2, ± ) form the final flavor state 
\up) at the detector 

I Up) = Up t+ \$ 2 ,tP2, + ) + ^-1^2,^2,-) • (12) 

Note that the initial and final flavor states are shown as red bullets in Figs. [I] and 0. 
I cpi, = 1,2) correspond to the orthogonal pairs of mass (energy) eigenstates lying 
on distinct elliptic axes depending upon the values of tfi (shown as blue bullets in Figs. [1] 
and @). 

To illustrate this point of view, let us carefully examine the expression for transition 
probability along the lines of Ref . [[| . The probability for flavor transition u a — > up is given 
by squaring the amplitude, 

V{u a ^up) = \A{u a ^up)\ 2 

= (U a \ <fl, + )( $2,^2, + I Vp)(up I 1? 2 , V?2, + V?l, + I ) 

+ ( U a I ^1, V9i, - )( -& 2 , <P2, ~ I Vp ){ Up I 7? 2 , t P2,-){A, ~\^a) 

+ [( u a I ^1, - )e iV - ( 7? 2 , ^2, - I ^ )( ^ I #2, <^2, + )e" iD+ ( 01, + I ) + C.C.] . 

(13) 

The two direct terms correspond to classical probability and two cross terms contain the 
effect of geometric and dynamical phases. In the limit when coherence is lost (for instance 
when the path length is much larger than the oscillation length), the total probability 
reduces to the classical probability. The quantum mechanical effect of oscillation essentially 
lies in the cross terms. 

At this juncture, we would like to introduce Panchar at nam's idea of closed loop quantum 
collapses and the related geometric interpretation of the cross terms appearing in the ex- 
pression of probability (Eq. (|1"3"D). Given any three rays | 21), | OS ) and | <£) in the projective 
Hilbert space for a two state system (such that the neighbouring rays are non-orthogonal) 
the phase of the complex number (2l|(£)((£|93)(93|2l)is given by Q/2 where Q is the solid 
angle subtended by the geodesic triangle at the center of the Poincare sphere. This excess 
phase is known as the Pancharatnam phase. Panchar at nam's phase reflects the curvature 
of ray space and is independent of any parameterization or slow variation. Thus it can also 
appear in situations where the Hamiltonian is constant in time i.e. for the case of vacuum or 
constant density matter in the case of neutrino oscillations. Furthermore, note the fact that 
Schrodinger evolution (possibly) interrupted by measurements can lead to Pancharatnam's 
phase. The crucial requirement for Pancharatnam's phase to be well-defined is to have 
cyclic projection of (at least three) states and the consecutive collapses should be between 
nonorthogonal rays. 



In Eq. (0) the cross term is related to the two path interferometer in energy space 



shown in Ref. M. The structure of the cross term (without the dynamical evolution) is a 
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series of cyclic quantum collapses with intermediate adiabatic evolutions given by | u a ) — > 
| tpx, + ) -» 1 1? 2) if 2, + ) -> I ^ ) -> 1 1?2> ^2, - } -> I Vi, - ) -> | ) which essentially 
maps a trajectory (as illustrated in Figs. [I] and 0) and subtends a solid angle of Q at the 
center of the Poincare sphere. The center of the Poincare sphere is a singular point which is 
the same as the degeneracy point corresponding to the null Hamiltonian. So, the appearance 
of geometric phase is intimately connected to the existence of this singular point. 

The interference term (without the dynamical phase) can be expressed as re* fi and 
therefore picks up a phase which will be B = Q/2 (half the solid angle) due to Pancharatnam's 
prescription. Note that even though the collapse process leads to a loss in probability (r / 1) 
the phase of the cross term remains unaltered whether we perform continuous evolution 
(unitary) or collapses (non-unitary) as long as it is cyclic. 

The Pancharatnam geometric phase is given by 

B = Alg[(u a | #i, (p u - )(l? 2 , <f2, - I Vp){vi3 I $2, V2, + )(#1> fl: + I V a )] (14) 

A few remarks concerning the states | v a ) and | up ) are in order. The states | u a ) 
and \ vp) can be any general flavor states. Pure flavor states (\v e ) and |^)) are very 
special (real) linear superposition of mass eigenstates and correspond to those produced and 
detected via weak interactions. If either \u a ) or \vp) are identical to any one of the four 
eigenstates involved before and after the adiabatic evolution, the interference term vanishes 
and then it is meaningless to talk about the phase of such a term, whether dynamical 
or geometric. It is easy to understand this from Pancharatnam's definition of parallelism 
between any two nonorthogonal states. 

For the case of appearance probability V [y a — > vp) one gets a geometric phase Q app = 
Q/2. For vp = u a , i.e. survival probability V [y a — > u a ) of the same flavor, it is easy to 
see that the path formed by cyclic collapses will again form a closed trajectory and subtend 
a non-zero solid angle at the origin. Therefore the interference term in survival probability 
will pick up a geometric phase which can be computed once we know the geometric phase 
in the case of appearance probability to be fl/2 by imposing unitarity which means that the 
interference terms should have a relative sign. Thus the cross term in survival probability will 
pick up a phase given by e lS = e*^-^/ 2 ). Hence both appearance and survival probabilities 
will pick up different geometric phases 

8 a PP = Q/2 and Q surv = n - fi/2 . (15) 

Next we discuss the distinct physical situations that can arise. In Figs. [I] and 0, we pictorially 
depict the geometric differences in cross terms between the CPC and CPV cases both for 
constant and varying density situations. The CPC case was discussed in Ref. D. In the 
present study with CPV phases, the density and CPV phase are two parameters which can 
either be constant or varying with distance, we list three possible cases : 

(A). Constant density case :- The pictorial depiction of the direction of collapse pro- 
cesses is shown in Fig. 0. In this case, we have four states (the two pair of mass and 
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Figure 1: Pictorial depiction of the cross terms (Eq. [Tj|) in case of appearance and survival probability 
for the CPC case. The direction of collapse processes connecting the four (six) neutrino states appearing in 
constant (varying) density case in the cross term are represented on the Poincare sphere by straight lines 
with arrows (see text). Note that the CPC mass states are labelled by | ± ) (i = 1,2). Cases (a) and (b) 
correspond to constant density situation while (c) and (d) stand for varying density situation. The states 
|i?2ji) are obtained after adiabatically evolving | i?i,±). Adiabatic evolution between the mass states is 
shown by a dotted line while dashed lines refer to collapse processes. 
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Qapp = n Qsurv = Q 




(c) (d) 

Figure 2: Direction of collapse processes connecting the four (six) neutrino states in absence (presence) 
of variation of density (and CPV phases) corresponding to the cross terms in probability (Eq. ([j~3])) are 
represented on the Poincare sphere by straight lines with arrows (see text). Note that the in presence of 
CPV, | &i,ipi,±)(i = 1,2) are used to denote the mass states. Cases (a) and (b) correspond to appearance 
and survival probability in a CPV (<p ^ 0) but constant density situation. Cases (c) and (d) stand for 
appearance and survival probability in presence of CPV and varying density case. Adiabatic evolution 
between the mass states is shown by a dotted line while dashed lines refer to collapse processes. 
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flavor states being orthogonal) that will lie on a great circle as shown in Fig. |2| (a) and 
(b). By appropriate gauge transformations the great circle can be made to coincide 
with the equatorial great circle. The bottomline is that we can use real form of the 
leptonic mixing matrix (Eq. (|9|)) and the phase will remain topological as in the CPC 
case . The Pancharatnam geometric phase is given by 

6 = Arg[±cos 2 (^/2)sin 2 (?9/2)] = 7t, . (16) 

Using Pancharatnam's prescription of quantum collapses, we have shown that in the 
two flavor context, CPV phase leaves no imprint of its existence in oscillation proba- 
bility for the case of constant density matter, hence providing a purely geometric view 
of this result. 

It is then natural to ask if CPV plays any relevant role in presence of adiabatic evolu- 
tion or alternatively one can ask if it is possible to gauge away the CPV phase which 
itself is evolving in time (or distance). We again use Pancharatnam's prescription to 
analyze this situation. 

(B). Density and CPV phase varying with distance :- In this case, we have the com- 
plex Hamiltonian (Eq. ([!])) and the eigenstates will not be real in general so the ray 
space is the full Poincare sphere S 2 (instead of a great circle S 1 ) in contrast to the 
case (A) listed above (see Fig. ^| (c) and (d)). There are four inner products appearing 
in the cross term, ( v a | ip l7 - ), ($ 2 , (p 2 , — \vp), (vp\ $2, f2, + ) and (#1, (pi, + \u a ) 
which are connected to the elements of the mixing matrix 0. We saw that by a 
suitable gauge transformation, the mixing matrix can be made real for two flavor case 
(Eq. (D). This implies that | u a ) and | $1, <£>x, ± ) can be chosen to be real and ly- 
ing on the equatorial great circle (x-z plane). However once this transformation is 
carried out, we do not have the freedom to choose the phases of the mass and flavor 
states any further. In general, the adiabatically evolved mass states |-$2,y?2,±) will 
be complex. Also, the final flavor state | vp ) which will be a linear superposition of 
states |-#2,y?2,±) with complex coefficients given by Eq. (0). The three states | vp) 
and | $2, ^2, ± ) will also lie on a different great circle and this is clearly illustrated 
in Fig. ^| (c) and (d). The fact that the two great circles (one containing the initial 
flavor and mass states and the other one containing the final flavor and mass states) 
do not coincide will lead to a net geometric phase (which differs from topological phase 
obtained for the CPC case 0). 

Let us compute explicitly the cross term for this case. Two of the inner products 
are real and positive, i.e., ( v a | $1, <fii, — ) = sin^i and ($1, <pi, + \ v a ) = cosi?i. 
But the other two are complex, given by ($2, — \ vp) = cos$2e** _Q+/3+7 - ) and 
( i/p 1 1?2) <P2, + ) = — sin -$2e - * 7 . This implies that in general the geometric phase picked 
up will be given by 

Qapp = n _( a _^ and fisurv = ^ _ Qapp = ^ _ ^ _ ^ 
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In the limit, (a — 0) — > 0, the phase becomes topological and that would actually 
correspond to a constant density (CPC or CPV) situation. 

Geometrically, we see that during the course of adiabatic evolution the mass state 
| $ 2 j <fi2, =t ) gets lifted from one plane containing \ $ 1 ,{p 1 ,±) and thus explore the 
full ray space instead of a great circle (with ( $2, y?2, ± | $1, <pi, ± ) being complex in 
general) . We have shown above that when presence of CPV induced by background is 
coupled with an adiabatic evolution with varying background density, the cross term 
of the probability (see Eq. (p~3|) ) will depend upon the extent of CPV or variation of 
ifii under adiabatic evolution (<^2 7^ <^i)- This leads to a geometric phase. Turning 
this argument around, we can conclude that CPV in the Hamiltonian can cause a 
deviation of e l& from ±1 in presence of an adiabatic evolution. 

(C). Density varying and CPV phase constant :- In this case, all the six states in- 
volved in the collapse process with intermediate adiabatic evolutions will be confined 
to a great circle and the phase will be topological in nature. This will be identical to 
the CPC situation. 

The key result of this article can be stated as follows. The topological phase obtained in 
Ref. for CPC situation was robust and did not depend up on any details of time evolution. 
We have shown that the presence of CPV in the Hamiltonian leads to a deviation of e 
from ±1 in presence of an adiabatic evolution iff the CPV phase also varies in space or in 
time. 

4 The cross terms and the action of CP 

Let us examine the form of cross terms appearing in general expression (Eq. (|IBD) for 
probability 

% = [( fa I #1, <Pl, ~ ( 2 , V*,-\vp)(vp\ 02, <^2, + )e~ W+ (A,<Pl, + \v a ) + C.C.] 

= 2rcos[6 + (XL -£>+)] 

= 2r [cos fi cos(XL - V+) - sin fi sin(£L - V + )} , (18) 

where r is the (real) amplitude of the cross term and r < 1. Geometrically we can see that 
the amplitude r will not change under CP. Inner product of two rays is a measure of distance 
between the two rays and this distance will not change under CP (reflection) operation. 

From the above expression (Eq. (|T%D), it is clear that for the special case when fi = ir or 
i.e. when the phase is topological in character, sinfi — > and we have only the first term 
(cos 6 term) contributing to the probability. However, when the phase becomes geometric, 
we will have both the terms contributing to the cross term. 
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Now we discuss the effect of CP transformation on the expression for transition prob- 
ability (Eq. flT3|)). CP transformation is an anti-unitary operation like the time-reversal 
operation and is given by tpi —> —ifi (i =1,2) 

Z CP = [{ Ua \ ti u -V?!, - }e W - (0 2 , ~<P2, ~ | Vf) )( V | 2 , -^ 2 , + )e"^ + ( #1, -Pl, + K ) + C.C.] 
= ^ re -i& e i(V--v + ) _|_ re jfi e -i(X>_-D + )j 

= 2rcos[(-fi) + (D_ -£>+)] 

= 2r[cosficos(£L -£>+) + sinfisin(£L -£>+)] . (19) 

CP transformation connects two different Hamiltonians or two different backgrounds that 
are CP partners, connected by (pi — > —ifi (i — 1,2). If we look at the form of cross term 
(Eq. fll8|)) and CP-transformed cross term (Eq. (|I9D), we note that the cross term and 
hence the transition probability will in general be different for the two CP partners. For 
the special case, when the phase is topological (sinfi — > 0) irrespective of the CP properties 
of the medium, the probability is invariant under CP. The cross term will no longer remain 
invariant if the presence of CPV is coupled with an adiabatic evolution as the phase fi 
becomes geometric. Pairwise the two Hamiltonians (connected by CP transformation) are 
distinguishable at the level of probability and CPV coupled with adiabatic evolution plays 
a crucial role for this to happen. 

We now discuss the implications of Wigner's theorem and ray space isometries in the 
present context. We can draw a connection between the above results and those rigorously 
obtained in Ref. |37j for any quantum mechanical system. It was pointed out that Wigner's 



theorem and the isometries of ray space can be most easily understood by using geometric 
ideas. The Pancharatnam phase fi itself is not an invariant under isometries of the ray 
space, but cosfi is. So under isometries of the ray space, fi can be preserved (unitary 
rotations) or reversed (anti-unitary reflections). Wigner 8 proved that any ray space isometry 
can be realized on the Hilbert space of quantum mechanics by either a unitary or anti- 
unitary transformation. The only measurable quantity is the transition probability which 
is defined as an overlap between rays and measures the distance between them. This result 
is captured in Eqs. (|18|) and (|I"9|). In two flavor case, when the geometric phase is restricted 
to be topological (CPC case and CPV in constant density situation), the cross terms in 
the probability depend only on cosfi and we can not distinguish between unitary (6 — > 6) 
and anti-unitary (fi — > —6) operations. But, in case the phase is geometric as in case of 
adiabatic evolution coupled with CPV, we have an extra term oc sin fi which will change sign 
under CP transformation if CP is violated. In our present work, we have given a geometric 
visualization of the operation of CP in case of two flavor neutrino oscillations, which is anti- 
unitary. Thus, for the two state system, the Poincare sphere description is quite rich and 
it encompasses the full set of possible operations, including unitary as well as anti-unitary 
operations. 



3 A complete algebraic proof of Wigner's theorem is given in Ref. 
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5 Conclusion 



We have shown that inclusion of CPV in the two flavor neutrino Hamiltonian can make 
the Pancharatnam phase geometric (and path-dependent) from topological, in case of adia- 
batic evolution of the mass eigenstates only when the CPV phases also vary with distance. 
The standard quantum mechanical treatment of two flavor neutrino oscillation contains 
Pancharatnam's phase which no longer remains topological when varying CPV phases and 
varying background density are taken into account. We have given a neat geometric illus- 
tration of the fact that the two flavor oscillation probability is sensitive to the presence of 
CPV phase in the Hamiltonian for the case adiabatically varying density profile. 

Concerning the observability of the geometric phase, we should note that the cross term 
in the probability should be accessible to measurements. We mentioned two physical settings 
in Sec. Q where the CPV phase could appear in the Hamiltonian. 

In the first case of ordinary matter and CPV induced by NSI in propagation, we should 
note that the extrinsic CPV phases are due to the e a o terms, which do not vary with distance. 
This would correspond to case (A) (if density is constant) or case (C) (if density is varying) 
described in Sec. |||. In both the cases, we have shown that the phase remains topological. 
A special case (R-parity violating supersymmetry) of non-standard interactions in ordinary 
matter was pointed out in Ref. |39| to show that in two flavor case, it is possible to have 
more than one essential parameters (electron, proton and neutron number densities) leading 
to non-zero Berry's jE| phase if the cyclic condition was also met in certain long baseline 



neutrino oscillation experiments. Note that there are other more exotic non-standard ways to 
get flavor off-diagonal and CPV terms leading to a modification of the dispersion relation in 



the neutrino Hamiltonian for example, violation of Lorentz invariance or CPT violation pT|| . 

The second case discussed in Sec. |^ was that of dense neutrino backgrounds where 
CPV can appear even when SI are taken into consideration and this was solely due to the 
fact that background neutrinos in mixed state can induce CPV terms in the Hamiltonian. 
Neutrinos traversing dense supernova cores in a region outside the neutrinosphere pass 
through varying electron and neutrino number densities 0. This could actually correspond 
to case B of Sec. if the CPV phase was also varying with distance where the Pancharatnam 
prescription predicts a geometric phase. In case of neutrinos coming from Supernovae it is 
impossible to detect the effects due to CPV induced geometric phases since the oscillations 
would get averaged out by the time the neutrinos reach the detectors on the Earth and the 
cross term will not be accessible to measurements. 

If in a hypothetical situation, it was possible to realise both varying CPV phases and 
varying density, then Pancharatnam's prescription predicts a geometric phase in the transi- 
tion probability different from n, 0. The problem then reduces to disentangling the geometric 
component from the total phase in the transition probability. The dynamical phase depends 
on the energy while the geometric part is independent of energy. So, by carrying out mea- 
surements at different energy values, we can extract the geometric part from the total phase 
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measured in oscillation experiments. 

To conclude it turns out that with only two flavors the robustness of the topological 
phase is very hard to destroy even in the presence of CPV as long as the density and/or the 
CPV phases are constant. For the two physical settings discussed above, the phase remains 
topological whether CP is violated or not. Thus, it calls for a study of the full three flavor 
situation where the physics of geometric phase can play a role [[12j due to the presence of 
intrinsic CPV phase in the leptonic mixing matrix. 
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